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Abstract. This paper investigates some properties of entropy solutions of hyperbolic conser- 
vation laws on a Riemannian manifold. First, we generalize the Total Variation Diminishing 
(TVD) property to manifolds, by deriving conditions on the flux of the conservation law and a 
given vector field ensuring that the total variation of the solution along the integral curves of 
the vector field is non-increasing in time. Our results are next specialized to the important case 
of a flow on the 2-sphere, and examples of flux are discussed. Second, we establish the conver- 
gence of the finite volume methods based on numerical flux-functions satisfying monotonicity 
properties. Our proof requires detailed estimates on the entropy dissipation, and extends to 
general manifolds an earlier proof by Cockburn, Coquel, and LeFloch in the Euclidian case. 



1. Introduction 

In this paper, following pQ, we investigate some properties of entropy solutions to the Cauchy 
problem for hyperbolic conservation laws on a manifold : 

(1.1) d t u + V g -f(u,-) = 0, u — u(t, x) G R, t>0, x€M, 

Here, M is a d-dimensional, smooth manifold endowed with a Riemannian metric g and, for 
each constant u E R, the map x i— ► f(u,x) is a smooth vector field on M, i.e. a section of the 
tangent bundle TM. The well-posedness theory for the hyperbolic conservation law was recently 
established in Ben- Artzi and LeFloch ,1 . As in Kruzkov's theory [5] which applies to conservation 
laws in the Euclidian setting M = R d , one is interested in weak solutions of (|1 . ip in the sense 
of distributions that are constrained by entropy inequalities. For instance, in the case that / is 
divergence-free, that is, 

(1.2) V g -f(u, •) = for every constant u, 
the entropy inequalities read 

(1.3) d t U(u) + V g -F(u,-)<0, 

for every convex U : K — > M, where the vector field x >— * F(u, x) is 

d u F{u, x) := d u U(u) d u f(u, x), lel,i£ M. 

Let us emphasize that the equation (|1.1|) is a geometric partial differential equation which 
naturally depends on the geometry of the manifold, only. In particular, all estimates derived on 
solutions to (|1.1[) should take a form that is completely independent of any particular system of 
local coordinates of M, although in practice, for the proofs it will often be convenient to introduce 
a particular chart to represent the manifold. Throughout, the convention of implicit summation 
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over repeated indices will be used. In local coordinates x — (x J )i<j<d, we will use the short-hand 
notation dj := djdxK Recall that the divergence operator arising in (jl.ip takes the form 

Vg-.f(u(t,x),x) := 2—^ dj(y/\g(x)\ f (u(t, x), x)) , 

where (gij) are the coordinates of the metric tensor g and \g\ := det(gij). Here, f J (u,x) are, for 
every u, the coordinates of the vector held f(u, •) at the point x. 

The present paper supplements the well-posedncss results established in the companion paper 
PQ and has two main objectives. 

First, we investigate some properties of the total variation of solutions of (jl.ip on a manifold. 
Precisely, we generalize the Total Variation Diminishing (TVD) property to manifolds, by deriving 
conditions on the flux / and a given vector held X ensuring that the total variation of the solution 
along the integral curves of X is non-increasing in time. Recall that TVD schemes for nonlinear 
hyperbolic problems play a central role in scientific computation, for instance in gas dynamics. 
Our diminishing total variation properties provide certain a priori estimates which, for instance, 
can be tested numerically and may help in designing robust schemes that are consistent with 
large-time asymptotics. Note that solutions of a conservation law on a manifold need not have a 
total variation that is bounded uniformly in time. The geometric effets may contribute to amplify 
the wave strengths, and the total variation of a solution may blow-up in the large as t — * oo. We 
also investigate here in some detail the case where the manifold is the sphere S 2 embedded in 
R 3 . Examples of hux (as vector fields in R 3 tangent to the sphere) are discussed from both the 
intrinsic and embedded standpoints. 

Second, we consider the numerical approximation of the entropy solutions of (jl.ip and we 
establish the convergence of the finite volume scheme on manifolds when the numerical flux- 
functions are depend monotonically upon their arguments. Our result is an extension to general 
manifolds of a theorem due to Cockburn, Coquel, and LeFloch [2] (see also [3]) in the Euclidian 
case. While the convergence of the finite volume scheme is easily established on a cartesian mesh, 
due to the invariance by translation of both the equation and the mesh, the convergence proof for 
non-cartesian meshes is more involved and, as was pointed out in [2], requires DiPerna's concept 
of measure- valued solutions [4]. To handle (jl.ip . we therefore must rely on Ben-Artzi-LeFloch's 
extension to manifolds [T] of DiPerna's theorem. 

An outline of this paper is as follows. In Section 2, we derive conditions on the flux and the 
vector field X ensuring that the function X{u) satisfies a conservation law whenever u is a solution 
of (jl.ip . This leads us to a diminishing total variation property for a class of flux and vector fields. 
In Section 3, we discuss the structure and general properties of conservation laws on the sphere 
S 2 . In Section 4, we apply the framework of Section 2 to the case of the sphere and derive total 
variation bounds. Next, Sections 5 and 6 are devoted to the statement and to the proof of the 
convergence of the finite volume scheme. 

2. Total variation diminishing estimates on a general manifold 

In the present section, we derive conditions on the flux / of the conservation law (jl.ip which 
ensure that the total variation (at least along certain vector fields) of the entropy solutions of (jl.ip 
is non-increasing in time. We will first state the main results (Proposition 12. ll and Theorem 12. 3p 
and then recall some elementary notions from differential geometry, before giving the proofs of 
the results. 

2.1. Statement of the estimate. Throughout, M is a smooth, <i-dimensional, Riemannian 
manifold [M,g), which has no boundary and need not be compact. The following result provides 
us with a key identity on directional derivatives of solutions of (jl.ip . 

Proposition 2.1. Let u : R + x M — > R be a smooth solution of the conservation law (jl.ip on 

M where, f = f(u,-) is a (smooth) vector field depending on the parameter u. Then, given any 
(smooth) vector field X on M the function w := X(u) : R + x M — > R satisfies the (linear) 
hyperbolic equation 

(2.1) d t w + V g -(wf u (u,-))=-g(V g u,(£xfu)(u,-))-X(V g -f)(u,-), 
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where (Cxf u ){ u , ■) denotes the Lie bracket of the vectors fields X, f u , that is 

Cxfu{u,-)~[fu(%-),X], ueR. 

The first term in the right-hand side of (|2.ip depends on the gradient of the unknown function 
u while the second term depends on the function itself. Observe that our formula is non-trivial 
even when M is the Euclidian space W 1 . In that case, the equation (]2.ip becomes 

d t w + V(wf u {u, x)) = - Vu • [X, f u ] (u, x)-X- V(V/) {u, x), 

where w — X ■ Vtt. When the flux is independent of x it is obvious that for any constant vector 
field the right-hand side of (|2.1|) vanishes and we can recover the total variation diminishing 
property for scalar conservation laws in M. d . 

We denote by TV the total variation functional for functions defined on the Riemannian man- 
ifold (M, g) . We introduce here a generalization, which extends to discontinuous functions the 
formula, valid for smooth u, 

TV x (u) := [ \X{u{t,x))\dv g . 
Jm 

Definition 2.2. 1. To any function u and vector field X defined on the manifold M one associates 
the quantity 

TV x (u):=sup [ uV g -{(f>X)dv g , 

where dv g is the volume element on the manifold and the supremum is taken over all smooth 
functions 4> '■ M —* R satisfying ||<^||l°° < 1. When TVx(u) < oo the function u is said to have 
bounded total variation along X. 

2. A flux f — /(u, •) on M is said to be divergence- free if 

V ff -/(u,-) = 0, ueR. 



It follows from Proposition 12.11 and the existence theory in [T] that : 

Theorem 2.3. (Total variation estimates.) Let u be an entropy solution of the conservation law 
(jl.ip on the manifold M where f is a smooth flux. 

1. The total variation of u along a vector field X satisfies 

TV x (u(t))<TV x (u(0))+ sup \{L x f u ){u,-)\ g I TV{u{,-))dt 

(2.2) (o,t)xAf Jo 

+ ||^(V s -/)(«,-)|Ui((0,t)xM)- 

2. Hence, when f is divergence-free and the Lie bracket of f u and X vanishes, 

C x fu{u, ■) = 0, ueR, 

then the solution has bounded total variation along X for all times t > 0, if this property holds at 
the time t = 0, and moreover, 

(2.3) TV x (u(t)) < TV x (u(0)), t > 0. 

At this juncture, it is important to recall from [1] that entropy solutions of (jl . 1|) have bounded 
variation (in the standard sense) at every time t, if this is true at the time t = 0. Actually, one 
has the general estimate 

TV(u(t)) < e Cl 4 TV(«(0)) + C 2 , 
where C\,C2 depend on the metric g and on derivatives (up to second order) of the flux /. In 
contrast, our result in Theorem 12.31 solely assumes that TV x (u(0)) is finite and concludes that 
TVx (u(t)) is finite for all times. The property (|2.3p should be useful when designing approximation 
schemes for (|1.1|) . and it is natural to require (|2.3|) to hold at every time-step of the discretization. 

Let us now recall some basic notions and notations that will be used throughout this paper. 
The differential of a function u : M — > R is the field of 1-forms defined by 

du(X) := X(u), 
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for every vector field X. It is convenient to choose a basis of the tangent space T X M at a point 
x E M together with its dual basis associated with the cotangent space T*M, so that we can 
introduce coordinates for vectors and covectors relative to these bases. As usual, to each vector 
{X 1 ) s T X M we associate its covector {g^ X 1 ) by lowering indices using the metric (implicit 
summation on repeated indices being used). Denoting (g^) the components of the inverse matrix 
associated with (gij), then, to each covector (rjj) G T*M we can also associate the vector (g^ r]j). 

Recall also that the gradient of a (smooth) function u : M — > M is the vector field V g u 
associated with the field of differential forms du. The divergence of a (smooth) vector field / on 
M is the function V s -/ : M -> R defined by 



/ uV g -fdv:=- du(f)dv g 
Jm Jm 



i m Jm 

for all smooth u : M — > M, where dv g is the volume element on M. 

Consider now a coordinate chart (x l ) together with the associated basis of vectors {d/dx l ) and 
covectors (dx l ). The differential of a function u is given by 

du = dx l . 
ox 1 

In view of the relations V*w := (V g u) 1 = g^ (du)j, we see that the components of V s u in the 
basis djdx 1 of T X M are 

du 
dx^ 



If the support of a function u is included in the domain of definition of the coordinate chart (x*), 
recalling that \g\ denotes the determinant of the metric tensor we can write 



J M d U{f) d Vg = fJ^r)^...d X « 



Jm ^\dx^ J 



so that 



The expression of V 9 • / follows also from computing the covariant derivatives of the vector 
field, as follows. Given a local chart (x l ) for M, the covariant derivative of the vector field / is 
the (1, l)-tensor field Vfc/- 7 defined by 

V k f : ihP* ■ YiJ 1 . 

where T ] kl are the Christoffel symbols given by 



In particular, we have 



and, since di\g\ = \g\g kl dig k i, 



T ii '■= \ { d kgu + dig k i - d l g k i)g 13 ■ 



it- - y j o k9ij 
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The divergence of a vector field / is the trace of the covariant derivative, i.e. 

v fl -/ := v k f k = d k f k + rt fe f = d k f + -^(vIsDf 



151 

d k f k + ^=(- VWldkf + d k (V\i\f k )) 



which is 

2.2. First approach. The proof of Proposition 12. II will follow from the following two technical 
lemmas. 

Lemma 2.4. For every smooth function u : M — ► M f/ie following identities holds (x £ M) : 

(2.5) V fl •(/(«,!)) = du(f u (u,x)) + (W g -f)(u,x) 

(2.6) £ x (^(w,a;)) = X(u)y u («,a;) + (£ x F)(m,j;) /or a// vector fields Y = Y(u,x). 

(2.7) X(h(u, x)) = X(u)h u (u, x) + X{K)[u, x) for all functions h. 
Proof. In local coordinates we have 

W g -(f(u,x)) = ^=d i (^\g~\f(u,x)) 

= (V fl -/)(u,a;) + du(f u (u,x)), 
which proves (|2.5p . To prove (|2.6p we note that 

(Cx(Y(u,x))) 1 = [X,Y{u,x)] i = X j d j {Y\u,x))-Y\u,x)d j X i 

= {x'djY* - Y 3 djX l ) (u, x) + x J 'y^(u, a:)a,-« 

= [X, Y]\u, x) + X^djuY^u, x). 
This shows (|2.6p . The proof of (|2.T|) is completely similar. □ 

Lemma 2.5. For any (smooth) function u : M — > H., vector field X , and flux f = f(u,-), the 
following identity holds : 

(2.8) X{V g ■ (f(u, x))) = V 9 • (X(u)f u {u, x)) + g{V g u, (Cxfu)(u, x)) + X(V fl •/)(«, 
Proo/. Using (f2~5]) and (|2~7l) we find 

X(V g (f(u, a;))) = A" ((V ff •/)(«, *)) + *(<?(V g u, (/>, a))) 

= X(V fl ./)(«,a:)+^(«)(V fl -/i.)(« > a;)+A-( S (V fl ti, (/«(«, a:))). 

On the other hand, 

V 9 • (X («)/>, x)) = A» V ff ■ (/„(«, a:)) + g(V g (X(u))J u (u, x)) 

= X{u){V g -f u ){u, x) + X{u)g(V g u, f au (u, x)) + g(V g (X(u)), f u (u, x)), 
using (|2.5p . Thus, we arrive at 
X(V g -(f(u,x))) -V g \X{u)f u {u,x)) = X{V g -f){u,x)+X(g{V g uJ u {u,x)) 

- x ( u ) g(V 9 u, fuu(u, x)) - g(V g (X(u)),f u (u, x)). 
Expressing the Leibnitz rule for the Lie derivative we obtain 

X(g(y g u, f u (u,x))) = {£xg)(V g u,f u (u,x)) + g(£xV g u, f u (u,x)) + g(\7 g u,£ x (fu(u,x))). 
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Therefore, we have 

X(V g ■ (/(«, x))) - V fl - (X(u)f u (u, x)) = X(V g -f)(u, x) + {C x g)(V g u, f u (u, x)) 

+ g(CxV g u - V g (X(u)), f u (u, x)) 
+ g(y g u,C x (f u (u, x)) - X(u)f uu (u,x))). 

Using (|2.6p . the last term above equals g(V g u, (£xfu)(u,x)). 

The proof of the lemma will be complete once we check that, for all vector fields X, Z and 
functions u, 

(2.9) g(C x V g u - V g (X(u)), Z) = -(C x g)(V g u, Z). 
To this end, in local coordinates we can write 

(C x V g u - \/ g (X{u))Y = X 1 d 3 V l u - V 3 udjX l - V\X 3 Oju) 

= X 3 (OX'a - V*0jlt) - V j ud j X i - V'.V 0,«. 

For the first term above we have 

X j (djV'u - V l d 3 u) = X 3 (dj(g ik )d k u + g ik (djd k u - d k d 3 u)) 
= X{g' k )g kl Vu 
= -X{g kl )g ik V l u. 

For the remaining two terms, a tedious but straightforward computation gives 
-VludjX* - V*X j d jU = -^ug il {{Cx9)ji - X(g kj )), 
where the components of C x g are computed from the formula 

X(g(Y, Z)) = (C x g)(Y, Z) + g{C x Y, Z) + g(Y, C X Z). 

Therefore, we find 

(C x V g u - W g (X(u))Y = -g k ^u{C x g) ]k , 

and, finally, 

g{C x V g u - V g (X(u)), Z) = -g m g k ^u(£ x g) jk Z m 

= -Vu{C x g) jk Z k 
= -(C x g)(V g u,Z), 

which proves (|2.9p . This completes the proof of Lemma 12.51 □ 
Proof of Proposition \2.1[ Applying the vector field X to the conservation law , we find 

d t x( u ) + x(y g -(f(u,x))) = 0. 

The first term is precisely d t w. Using Lemma l2.5l to evaluate the second term, we find the desired 
identity (|2.1[) . This completes the proof. □ 

2.3. Second approach. We provide here a second approach to Proposition 12. 11 To any vector 
field X we associate the (local) one-parameter group of diffeomorphisms <p : M — » M (defined for 
all sufficiently small s) from the integral curves of the vector field X, 

(2.10) ^(x)=X i ( l p s (x)), <po(x)=x. 

as 

Given a function u : M — > R we may consider the composite function u o tp s . 

Lemma 2.6. Let X be a smooth vector field on M and (<p s ) its associated one-parameter group 
of diffeomorphisms. Then, for any smooth function u : M — ► M it holds 

V g -(f(uo(p s ,-)) 

[2 ' " ! = V g -(f(u, •)) o - s (x(V B 'f)(u, •) + g(V g u, (Cxfu)(u, •))) + 0(s 2 ). 



HYPERBOLIC CONSERVATION LAWS ON MANIFOLDS 



7 



Proof. We will first check that for any vector field Y, 

(2.12) Yi(ip a (x)) - d^Xix) = s£ x Y j + 0(s 2 ). 

First, it is obvious that 

Y ] a (p s = Y j + s X k d k Y 3 + 0{s 2 ). 

Second, observe that for all i,j 

^-d l ip 3 s = d^ifl = di(X j o tp s ) = (d k X 3 ) o ip s d l ip\. 
as as 

Also, since (fio( x ) = x, we have 8iyjl = 5{ and so ^9^1 _ = diXK Therefore, Taylor 
expanding yields us 

ditpi = 5{ +sdiX j + 0{s 2 ). 

In consequence we find 

Y j (<p s (x)) - d l iflY l {x) = sX i d i Y j -sY l d l X J + 0(s 2 ) 
= sL x Y3 +C(s 2 ), 

which is precisely (j2.12|) . 

Setting v := u o <p s , we start from (|2.5[) which, for the function v, reads 

(2-13) V g - (/(«,-)) = (V g -f)(v,-)+dv(f u (v,-)). 

Let us compute the last term of the above identity. In local coordinates, we can write 

[dv)i = div = di(u o ip s ) 

= {(dju) o <p„) dupi, 

and therefore 

dv{f u (v, •)) = ((dju) o <p s ) d.ipi fl(v, •) 

= ((djU) o ip s ) (diifl fl(v, •) - fi(u, ■) o ip s ^ + (djufi(u, •)) o ip s 
= ((dju) o ip s ) (di<pl f l u (v, •) - fi(u, •) o ip s ^j + du(f u (u, •)) o (p s . 

Taking (|2.12p into account and noting that 

((d jU ) o <p.) £ x fi(v, •) = d jU £ x fi(u, •) + 0{s) 

= g(V g u,Cxfu(u,-)) + 0{s), 

we arrive at 

dv{f u (v,-)) 

= -s (d jU ) o <p 3 (C x fi)(v, •) + du(f u (u, •)) o <p s + 0(s 2 ). 

= -sg(V g u, (CxfuKu,-)) +du(f u (u,-))o Vs + (D(s 2 ) 

= -sg(V g u, (£xfu)(u, •)) + V s •(/(«, ■)) o <^ s - (V fl -/)(«, ■) ° + ^(s 2 ), 

where we have also used (|2.5[) with the function u. 

Returning now to (|2. 13|) . using the identity for dv(f u (v, •)), and re-ordering the terms, we find 

V 9 -(/(«, •)) = V,.(/(u, •)) o <fs + (V 3 •/)(«, •) - (V g -f)(u, •) o Vs 

-sg(V g u,(£ x fu)(u,-))+0(s 2 ). 

Finally, observing that 

(V g -f)(v,x) - (V g -f)(v,<p s (x)) = - s ±(V g .f)(v,<p s (x))\^ + O(s 2 ) 

= -sX k d k (V g -f)(u,-) + 0(s 2 ) 
= -sX(W g -f)(u,-) + 0(s 2 ) 
concludes the derivation of (|2.11[) . □ 
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Proof of Proposition We provide here a second proof of this theorem. If u is a smooth solution 
to (jl.ip then, observing that dtv = (dtu) o ip s , we can rewrite in the form 

d t v + V g ■ (f(u, x)) oip s = Q. 

So, using (|2.1ip we get 

(2.14) d t v + V g -(f(v, •)) - -a (x(V g -f)(v, ■) + g(\7 g u, (C x f u )(v, •))) + 0(s 2 ). 
Now, recalling the conservation law (jl.ip and the equation (|2.14p . we obtain 

dt (kv - «)) + V 9 • (i (/(«, s) - f(u, x))) 
= -X(V g -f)(v, x) - g(V g u, (C x f u )(v, x)) + (D(s), 



(2.15) 
in which 



and 



lim — (v — u) = X(u) 
s^Q s 



lim ~(f(v, ■) - f(u, ■)) = lim ^ (/(«, •)) 



= /„(«, •)X i 9 i u = X(u)/„(u,-). 
So, letting s tend to zero in (|2. 15[) yields the desired identity (|2.ip . □ 



3. Conservation laws on the sphere S 2 

Following [1] we include here a discussion of basic properties of the conservation law (jl.ip in 
the case (so important in the geophysical applications) that the manifold M is the sphere S 2 . We 
discuss here the choice and properties of the flux in (jl.ip and we emphasize differences between 
the intrinsic and the embedded representations of the sphere. Certain properties will be closely 
related to the fact that the sphere S 2 can be embedded isomctrically in the Euclidian space R 3 . 

3.1. Spherical coordinates on S 2 . One convenient approach consists in parameterizing the 
smooth manifold S 2 punctured at the North and South poles, by using a single chart in spherical 
coordinates (ip,8), where ip G [0, 27r] denotes the longitude and 8 S (0,7r) the latitude. Set 

O := T 1 x (0,tt), 

where T 1 = [0, 2ir] with the identification of the points ip — and 2tt. By definition, a smooth 
function defined on Q is extendable to a smooth periodic function for all tp € M. Indeed, in all 
what follows, we always consider smooth functions defined on S 2 represented as functions of (ip, 8) 
outside the poles and extended to the poles by continuity. We imbed S 2 in M 3 according to 

x = x(ip, 8) = (cos ip sin 8, sin ip sin 8, cos 0) . 

The Riemannian metric tensor of the punctured sphere, when expressed in spherical coordinates, 
read 

I \ /sin 2 8 0\ 
{9ii) = [ o 

Recall that the inverse matrix of g = (gij) is denoted by (g 1 -'). 

The components of the gradient V ' g u of a function u : S 2 — > M in the basis (d/dip, d/dd) of 
T X S 2 are 

According to (|2.4p , the divergence operator on the sphere is given by 
(3.1) ^.f:^ 9 (l(r^8) + §- d (f^8)), 
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for every smooth vector field / on S 2 , where f v and f e are the coordinates of the vector / in the 
basis (d/d<p,d/d9). Similarly, the Laplace-Beltrami operator is defined for all u : S 2 — > R by 

1 ( d t du \ d /du 



( 3 - 2 ) A g u := V s • (V g u) = -— — - +7jhr sin( 

sin \ ay V sin <9y / 90 V 96* 

We now turn to the following conservation law on the sphere : 
(3.3) 9 t tt + V fl -/(u,-) = 0, (i,i)el+xS 2 , 

where the flux / : R x S 2 — > T5 2 is a smooth vector field. In view of (|3.1|) . when spherical 
coordinates are used the conservation law (13.31) takes the form 



(3.4) a t (u sintf) + y,0) S in0) + — (f e (u,tp,0) sintf) = 0. 

ay of) 

A periodic boundary condition is imposed in the variable (p £ [0, 2tt], that is, 

(3.5) u(t,O,d)=u(t,2n,0), 9e[0,n],t>0. 

On the other hand, no boundary condition may be imposed along (9 = and 9 = n. Observe that 
the equation (|3.4p is singular (sin 9 = 0) precisely at the poles. 

If u is a smooth solution of the equation (|3.4p , a natural question is whether u satisfies some 
additional conservation laws. We have already noticed that the Total Variation Diminishing 
property is available for solutions to conservation laws in M. d when the flux-function is independent 
of the spatial variable. We have seen also that, when writing the conservation law in local 
coordinates, the flux in general depends explicitly on the space variables. It would be meaningless 
to require that the flux is "independent of x" . A natural approach is to restrict attention to flux 
that are divergence-free, which in spherical coordinates is expressed as 

(3.6) 5- (/*(«, sm9) + —(f e (u,ip,9) sin 9) = 0. 

Recall that an entropy / entropy flux pair for the conservation law (|3.1|) is a pair (U, F) in which 
the function U : R — * R is arbitrary and F(u, x) = J" U'(u) J^f(u, x) du. 
For the sake of completeness we recall : 

Proposition 3.1. (Intrinsic description.) Consider a conservation law (II. ip on the sphere S 2 
with a divergence-free vector field f = f{u,x), and let (U,F) = (U(u),F(u,x)) be any entropy 
/ entropy flux pair. Then, smooth solutions u = u{t,x) of (|3.4[) - (|3.5[) satisfy the additional 
conservation law 

(3.7) d t U(u) + V g -F(u,-) = 0. 
In particular 

(3.8) j t jju(u(t,-))dv g = 0, t>0. 

Of course, for weak solutions, (|3.7p - (|3.8p will hold as inequalities, only. In spherical coordinates, 
equation (|3.7p reads 

d t (U(u) sm9) + —(F^(u,^9) sinO) + — (F e (u, <p, 6) sm9) = 0. 
dip 09 

Proof. Following [T] we multiply (|l.ip by U'(u) and get 

d t (U(u)) + U'(u)V g -(f(u(t, a), as)) = 0. 

Using (|2.5p with F instead of /, we find 

V g -(F(u(t,x),x)) = (y g -F){u,x)+du(F u {u,x)) 

U'{v)(V g -f a ){v,x)dv + U'(u)du(f u (u,x)). 

Again by (|2.5p . since / is divergence- free, the first term above vanishes and 

U'(u)du{f u {u,x)) = U'{u)V g -(f{u,x)), 
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SO that 

V g -(F(u(t,x),x)) = U'(u)V g -(f(u,x)), 

which proves (|3.7p . 

Furthermore, since S 2 has no boundary, the integral of the divergence of any vector field 
vanishes, thus yielding ()3.8|1 . 

It is illustrative to also derive the equation (|3.7[) using spherical coordinates. To simplify the 
notation we set 

d 

h l = f l sind, i = <p,6; d a = ——, a = <p,9. 

da 

We denote here by h % u , h^, h\ respectively the derivatives of h l with respect to the first, the second, 
and the third variables. Then, the equation (|3.7[) becomes 

d t (u sm9) + d lp (hv(u,<p,9)) + d 9 (h e (u,tp,8)) =0. 

Multiplying it by U'(u), we get 

= d t {U{u) sin 9) + U' (u^h*) + U'{u)d e (h 9 ) 

= d t (U(u) sm6) + U'{u){h% + Kdtpu) + U'(u)(h e 6 + h e u d e u), 

or equivalently, with H := Fsin6>, so that H l u — F* sin#, 

d t (U(u) sm9) + H*d v u + H e u d e u + U'{u){h% + h e g ) 

= d t (U(u) sine) + d v H v + d e H e + Q = 0, 

where 

Q = U\u){h% + hl)-H*-H 6 6 . 

Suppose now that the condition (|3.6p is satisfied, that is, + h\ = for all u e R. Then, we 
have 

H$+H 9 g = U'(v)(hl v + h%)dv = U'(v)(h% + h$) u dv = 0, 
Jo Jo 
and therefore Q = 0, which completes the derivation of (|3.7p in spherical coordinates. □ 

3.2. Embedded representation of S 2 . In the embedded approach, the sphere S 2 is regarded 
as a submanifold of the Euclidian space R 3 , via some isometric embedding h : S 2 — » R 3 . The 
canonical Riemannian metric g on S 2 is induced by the canonical metric g on the Euclidian space 
R 3 , by g — h*g, where, by definition, the pull-back h* maps the 2-tensor field g on R 3 to the 
2-tensor field g on S 2 , as follows: 

(3.9) (h*g)^,r,):=g(h^,Kr,), ^rjeTS 2 . 
The tangent vector h+t; is the push-forward of £ £ TS 2 into R 3 : 

(3.10) (h*0(f) '■= £(/ ° h ) for ever y function / : R 3 -> R. 

From now on we fix some local coordinates (x 4 )i=i,2 on S 2 (for instance, the spherical coordi- 
nates described above) as well as some orthonormal coordinates (x a ) a =i.2,3 on R 3 . In particular, 
we have = 8 a p (the Kronecker symbol) . 

We have 

(3-n) %w = §Je, 

so that we may compute the expression of the metric tensor g in these local coordinates: 

dx l dxi ' 
where x a := g a fix^ . From this, it follows that 
, n -, n X dx a dx a 



HYPERBOLIC CONSERVATION LAWS ON MANIFOLDS 



11 



In other words, we have {gij) = J* J, where J = (dx^/dx 3 ) is the Jacobian matrix of the 
embedding h. 

In what follows, all vectors are considered to be in R 3 , unless otherwise stated. A vector tangent 
to the submanifold S 2 of R 3 is regarded as an element of R 3 through the identification 

T X S 2 C T X R 3 ~ R 3 . 

Consider now the spherical coordinates introduced earlier. The outward unit vector to S 2 in IR 3 
at the point x = (ip, 9) is given by 

h((p,6) = n(x) := (sin 9 cost/;, sin sin p, cos 9). 

At any given point x — (ip, 9), a normalized basis of the tangent plane to S 2 in R 3 (now viewed 
as a two-dimensional subspace of R 3 ) is given by the two vectors 

/ \ 1 9n , , N 
t H a; ) := = {-smf,cosf,0), 

sm 9 dip 

t 2 (x) := — — = (—cos 9 cosf, — cos 9 sin if, sin 9) . 

Given a vector £ in TS 2 with (intrinsic) spherical coordinates (£ v ,£ ), we associate to it the 
vector £ G R 3 defined by (|3.10f) . The components of £ in the basis (ti, t 2 , n) of R 3 are denoted by 
(f )k,<3. The vector £ is tangent to the sphere: using (|3.11|) to compute its coordinates, we find 

f=|V + |^ = sin^t 1 -^t 2 , 

so that in spherical coordinates the intrinsic and embedded coordinates of a tangent vector are 
related by 

(3.13) I 1 = smO^, i 2 = -£, e , f 3 = 0. 

Observe that the normal vector is well-defined for all 9 (and p). However, the North (9 = 0) 
and South {9 = tt) poles are singularities for the tangent vectors in spherical coordinates. Namely, 
the particular value of if is irrelevant when 9 = or 7r; however, the tangent vectors to S 2 do 
depend upon if when 9 = or 7r, so that a basis of the tangent space is not uniquely specified at 
the poles. This is consistent with the well-known fact that any globally smooth, tangent vector 
field on the sphere must vanish at least once. 

In (embedded) spherical coordinates, the flux / : R x il — > R 3 has the general form 

f(u, ip, 9) f\u, ip, 9) ti{<p, 9) + f 2 (u, p, 9) t 2 (p, 9). 
Note that / is tangent to the sphere. In view of (|3.ip . (|3. 13|) . we have 

p-=wa.0F, f 2 = -f, 
and the divergence operator takes the form 

v '-' ! -5b(£-»C*" 

In terms of the unknown function u = u(t,p,9), (|3 . 3|) is a conservation law with non-constant 
coefficients posed in the rectangular domain (t, p, 9) € R + x Q; 

(3.14) d t (u sin9) + —{f\ u ,p,9)) - —(f 2 ( u ,p,9) sinfl) = 0. 

For the sake of clarity we will always state the final results in the intrinsic notation f v , f , but 
in the course of the calculation it will be convenient to use the embedded notation f 1 , f 2 , which 
has the advantage to be based on an orthonormal basis of R 3 . 

Let us now define the tubular neighborhood T C R 3 of S* 2 consisting of all spheres with radius 
1/2 < r < 3/2, so that S 2 corresponds to r = 1. Then, any smooth tangent vector field f(u,x), 
x G S 2 , may be represented at each point of the sphere by 

(3.15) f(u,x) = n(x) x $(u,n(x)), 

where $(it, x), x G T is a smooth vector field in T C R 3 . We can thus write Q(u, x) — n(x)). 
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Lemma 3.2. With the notation (|3 . 1 5[) the equation (|3.14[) takes the form 

(3.16) d t (u sinfl) + A($( U ,„ M ) . _ ^($(„,n(^(9)) ~) = 0, 

Proof. Setting 

(sin cos — sin <p — cos 6 cos 
sin sin <p cos y> — cos 9 sin 
cos 9 sin ( 

we have 

n = A ^oj , ti = A fij , t 2 = 4 , 

and A -1 = A*. Since det(A) = 1 the basis (ti,t2,n) is positively oriented and 

n x ti =t2, n x t2 = — ti. 

Using coordinates, we have 

&(u,n(cp,9)) ^^(u.n^^-ti^e) 

and 

f(u, ip, 9) = n(<p, 9) x n(<p, 9)) tifo 9) + $ 2 ( u , n(ip, 9)) t 2 (cp, 0)) 

= ¥(u,n(<p,6))t 2 (<p,6) - $ 2 ( U ,n(^0))ti(^,0), 
so that f 1 = -$ 2 , / 2 = I 1 . But 

sm6\ dcpJ' V 36 J' 

and so 

J V 8fi/' ^ sinfll 
or equivalently using the intrinsic notation 



dip J ' 
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sin6» V ' dfl/' ^ sin0V d<^' 
which, in view of (|3.14p , leads to (|3.16[) . □ 

Observe that, for smooth solutions, (|3.16[) takes the equivalent form 
(3.17) d t (u smff) + nfe *)).(__-__) = 0. 

We have the analogue of Proposition 13. II 







Proposition 3.3. (Embedded description.) Suppose that the vector field $ satisfies the divergence- 
free condition: for all u S M and (ip, 9) G O, 

(3.18) l^m^O)) • g M ) - |(*fcn&,,*)) • |( 

or, equivalently, 

< 3 " 19 > ^M-^-^ = °- 

Le£ [/ : K — » M and f : R x 5 2 ^ I 3 &e smooth maps such that 

d 

— < 

du 

Then every smooth solution u of (|3 . 1 6[) satisfies 

(3.21) sin0) +JL(iS(u,n(<p,0)) ■ - ^(*(«,n(M)) • ^) = 0. 



(3.20) V(u,n((p,0))= I U'(u)-^-^(u,n{ip,9))du. 
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Proof. Multiplying the equation (|3.16[) by U'(u) we obtain 

d t (U(u) sin 6) - d u V(u, n(<p, 8)) -d 6 u-^ + d u V(u, n(<p, 9)) -d v u-^ = 0, 

thus 

d t (U(u) sin 8) - d e (*(«, nfo 9))~^+d v {^{u, n{<p, 6)) • ^) + Q = 0, 

where 



dn dn 



U'(v)$ ue dv™- J U'(v)$ uv dv^ 



by the condition 1(5715)) , □ 

Remark 3.4. TTie constraint (|3.18|) can aZso oe written in the form 

dn\ dn / <9n\ <9n 



£)$ £)$ = o. 

In particular, if the map $ coincides with the gradient of some potential h, that is, $ = £>/i /or 
/i : T — > R, t/ien i/ie divergence- free condition (|3.18p is automatically satisfied. 

4. Total variation diminishing estimates on the sphere S 2 

We are now in a position to discuss the total variation of solutions to scalar conservation laws 
on the sphere. We are going to investigate the conditions discovered in Theorem 12.31 After 
comparing the intrinsic and embedded description it turned out that the second approach led to 
simpler (but still not so straightforward) and more intuitive formulas. In view of Section [3l a 
vector field tangent to the embedded sphere S 2 can be expressed as the cross product between 
the exterior normal n(x) and some other vector field which can always be chosen to be tangent 
to the sphere. So, given arbitrary flux / and vector field X let us introduce (smooth) maps 
$ : R — > T — > M 3 and * : T — > M 3 (which are not unique) such that 

f(u,<p,0)=n(<p,0) x$(u,n(<p,0)) 
N '" X(<p,9) = n(<p,9)x*(n(<p,0)) 

(where T C R 3 is the tubular neighborhood of S 2 defined in Section [3]). It should be kept in mind 
that all the vectors under consideration here belong to R 3 . 

Having in mind the class of divergence-free flux, we introduce : 

Definition 4.1. The flux f and vector field X are said to satisfy the gradient condition if there 
exist functions a : R x R 3 — > R and b : R 3 — > R depending solely on n = x/\x\ and are smooth in 
a neighborhood of the 2-sphere such that 

(4.2) $(u,x)=D x a(u,x/\z\), *(ar) = D x b(x/\x\), x e R 3 . 

From the gradient condition we deduce that, for instance, 

$ = D x a = D n a ■ D x (x/\x\). 

Clearly, D x (x/\x\)-n vanishes identically, since it is the derivative of the map x/\x\ in the direction 
of the vector n. Therefore, 

$ • n = * • n = 0, 

that is, the vectors <& and ^ are tangent to S 2 . Observe also that as in Remark l3.61 the divergence- 
free condition is automatically satisfied for fields satisfying the gradient condition. 
Here is the main result of this section. 
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Theorem 4.2. Let f and X be a flux and a vector field defined on the sphere S 2 respectively, 
satisfying the gradient condition (|4.ip ~ (|4.2p for some maps $,'1'. Then, the flux f is divergence- 
free and one has 

(4.3) £xfu(% -) = if and only if f u (u, ■) = C(u, -)X, u e R, 

for some scalar-valued flux-function C — C(u, •) satisfying X(C(u, •)) = 0. 

When (|4.3p holds, and when (V ff -f)(u,-) — 0, the conservation law (|3 . 14[) decouples into a 
one-parameter family of independent, one- dimensional conservation laws. 

Hence, in the case of flux which are "gradients" , we are able to identify the class of flux 
satisfying the condition of Section For each u, the vectors f u (u,-) and X are parallel and 
C(u, •) is constant along the integral curves of the vector field X. The curved geometry of the 
sphere restrict the class of flux enjoying the TVD property. This should be compared with the 
case of equations in the Euclidian space M. d where no such condition is implied on the flux. 

Throughout the rest of this section, a value u 6 R is fixed and so we omit it. Recall from 
Section [3] that the components (J 1 ,/ 2 ) of / in the basis (ti,t2) of the tangent plane to S 2 at a 
point x(ip, 9) read : 

f 1 = -$(n) • n , ~f = -r^$(n) • n v , 

sin a 

where we have used subscripts to indicate differentiation. Using the equation (|3 . 1 3[) which connects 
these components with the intrinsic components (/jf , /„) of f u , we find 

F = -^*(n) • n e 
sinC 



f d = —^Hn)-n v . 

SHlt) 



A similar formula holds for X 



X* = --^(n) • n e 

sm9 

X 9 = -—J>{n).n v . 
sin 



Observe now the following elementary identities: 
(4.4) 





— sin' 




cos 9 


n <p9 — 




sin 9 


nee = 


n. 



We are now in a position to reformulate condition Cxfu = 0. 

Lemma 4.3. The condition Cxfu = is equivalent to the requirement that = ^(n) satisfies 
the following system of partial differential equations with unknown X : 

($„ • ne)(^ v ■ n e ) - ($« • n v )(* e ■ n e ) = -(* • n v )({<f> u9 ■ n e ) - ($„ • n)) 

-(*-n)($»-n v ,) + (*-n 6 i)($ 
^ • n e )(^ v ■ n v ) - ($„ • n ¥ ,)(* e • n v ) = (* • n e )((<f> uip ■ n v ) - sin 2 0($„ ■ n)) 



+ sin 2 0(* ■ n)($ u • n e ) - (* • n ¥ ,)($„ e ■ n v ). 
Proof. The condition on the first component, (Cxfu) (p , of the vector Cxfu reads 

{X*d v + X 6 d B )ft - (ffd v + f u d e )X* = 0, 

that is 

— (* ■ n e )d v {<f>u ■ n e ) - (* • (-— ($„ • ne)) 

sinf/ Vsinff / 

-^-($ ■ ne)5 v (* • ne) + (4 • n v )d 6 (^—{^ ■ n e )) = 0. 
sin a \ sm t/ / 
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This equation takes the form 

(* • n e )($ ntp • n g ) - (* • n«p)($ U(0 • tie) + ' n v)( $ 

COS $ 

- ($ u • n e )(* v • n e ) + ($„ ■ n v )(* e • n e ) - (*„ • n ¥ ,)(4- • n e ) 

sinf 

+ (* ■ n e )($„ ■ n v g) - (* ■ n v )($„ • n ee ) 

- ($ u • n e )(* • n ev ) + ($„ • nip )(* • n ee ) = 0. 

In view of (|4.4[) . all terms above which contain the factor (cos 9/ sin 8) cancel out. Using that 
rigg = — n and re- arranging the terms leads to the first equation. 
The second equation is obtained similarly, by developing 

(C x f u ) e = (X*d v + X e d e )f u - (/„% + f u d e )X 9 

and using (|4~4]) . □ 

We now show: 

Lemma 4.4. A point x £ S 2 being fixed, there exist non-trivial flux f and vector fields X 
satisfying the gradient condition, and possibly defined in a small neighborhood ofx only, such that 
C-xfu = 0. These vector fields are characterized explicitly by the condition 

Xxf u = Cn 

for some scalar C = C(u) (independent of x). Among these solutions, the only flux f and vector 
field X that are globally defined and smooth on the whole sphere S 2 are given by 

f u (u,-) = C(u,-)X, 

where C = C(, •) is a (scalar-valued) function such that X(C(u, •)) = 0. 

Vector fields that are defined locally on S 2 may still be useful to control the total variation 
in a domain of dependence, but they do not provide information on the large-time behavior of 
solutions. 

Proof. First, let us show that for i = cp,6 we have a; = 4> • and bi — \& • m (subscripts indicate 
differentiation) . A straightforward calculation shows that 

x \ ( Id — n ( 

D x a = D n a ■ D x -j— r = D a a 

In particular, restricting the formula to the sphere we find D x a = D n a- (Id — n®n), and therefore 

di = D n a ■ rii = D n a ■ (Id — n <g> n) + D n a • (n (g) n) • nj = D x a ■ rii = <J> • n, , 

because (n <S> n) • m = 0. In consequence, we now simply write Da, Db. Recalling that by 
assumption we have Da ■ n = Db ■ n = 0, we can rewrite (|4.5p in terms of b: 

a ue (d v Db ■ tig) - a u Jd e Db ■ tig) = -bJd e Da u ■ tig) + bg(d v Da u ■ tig), 

(4.6) 

a u e(d v Db ■ ti v ) - a utp (d e Db • n^) = b g (d v Da u ■ n 9 ) - b ip (d g Da u ■ ti v ). 
This is now a system in the scalar variable b. Furthermore, using (|4.4|) we find 
a u 98 = dg(Da u ■ n g ) = d s Da u ■ n g - Da u ■ n = d e Da u ■ tig, 
CLiuptp = d v Da u ■ ti v + Da u ■ n w = d v Da u ■ n v — sin 9 cos 9 a u g , 

cos 9 cos 9 

due<p = o v Da u ■ tig + Da u ■ tt v g = o v Da u ■ tig + — — -Da u ■ tt v = o v Da u ■ tig + — — -a v 



sin f sin f 



so that 



d g Da u ■ tig = 
d v Da u ■ 11^ = a U(p(p + sin 9 cos 9 a. b 

on _ cos9 
o^JJa u ■ tig — a u g v : — -a uv . 
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Similarly, we obtain 

dgDb ■ n e = bg 



d v Db ■ n v = b w + sin 9 cos 9 bg 



d v Db ■ n s = bg v r-jiv 



cos6> 
sm9 

and, therefore, the first equation of the system (|4.6|) becomes 

cos 9 

Om a u gg + b v g a u g — a uip bgg — bg a uip g — — — - 

sm u 

This is equivalent to saying 



{a u g b v — bg a u(p ) g — ^ (a u g b v — bg a uip ) , 
which, after integration, leads to 

a u e b v - bg a uip = sin 9 e G ^ 

for some function G = G(ip). 

As for the second equation in (|4.6[) . we have 



a u e + b v a u g v - bg a uipip - a uip b v g + sin 9 cos 9 a u gbg - sin 9 cos 9 bga u g 

- 22ll( h -h \ 

so that 

(a u e b v — bg a uip ) = 

and a u g b v — bg a ulp is a function of 9 alone. Therefore, G = G(ip) is a constant and we find 

a u e b v — bg a uv = C sm9 

for some scalar C = C(u). 

The condition above can be written in terms of the original components of f u and X, yielding 

C 

Ju Ju sin0 

and, for the embedded components, 

'fix 1 - fix 2 = a 

Since X and /„ are both tangent to the sphere, the last equation means that there exists a 
constant C such that 

X x f u = Cn 

for all x £ S 2 , so 

\X x f u \ = C, xe S 2 . 

This formula describes vector fields that satisfy the constraint but may not be globally smooth. 

Imposing now that X is smooth, we deduce that it must vanish at least once on the sphere, as 
must the flux /„. By continuity, we necessarily have C = so that f u and X are colinear. This 
means that there exists a function C(u, x) such that f u — CX for all x S S 2 . Finally, substituting 
f u for C X in Cxfu = yields X(C) = 0. This completes the proof of Lemma T4.4I □ 

Proof of Theorem \4.£\ The first claim of the theorem is an immediate consequence of Lemma 14.41 
The equations are written along any integral curve of the vector field X. To see this, fix a 
point x € S 2 and introduce a parametrization s i— » p(s) of the (unique, integral) curve passing 
through x and such that the velocity vector -^<p(s) equals the vector X for every s. Then define 
the function 

u(t,s) :=u(t,tp(s)), t > 0, set. 
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Note that since X{C) = we have C(u(t,s),tp(s)) =: C{u) for all s. Therefore, using (|4.3|) and 
31). we deduce that 



d s C(u(t,s)) = C u (u(t,s))X(^(s))(u) = du(C u (u(t,<p(a)))X(<p(a))) 
= V 9 •(/(«(£, </j(s)),^(s))) = -dtu(t,ip(s)) = -d t u(t,s), 

thus 

9 t u(t, s) + S,C(u(t, *)) = 0. 

Note that this is a one- dimensional conservation law in which the flux C does not explicitly 
depend on the spatial coordinate. By Theorem l2.31 these conservation laws admit total variation 
bounds of the form (|2.3[) . provided that / is divergence- free. □ 

5. Finite volume method 

The rest of this paper is devoted to establishing the convergence of the finite volume method 
on a d-dimensional Riemannian manifold (M, g) . We will construct approximate solutions to 
the Cauchy problem associated with the conservation law and the initial condition [uq G 

L 1 (M)nL°°(M)) 

(5.1) u(0,x) = u (x), x G M, 

using a finite volume method based on monotone numerical flux-functions, and we will establish 
that such approximations converge strongly to the unique entropy solution of the problem. We 
need not assume that the flux / is divergence free (see (|1.2[1 ), and on the map V s -/ we require 
the following growth condition 

(5.2) \V g -f(u,x)\<C 1 + C 2 \u\, (u, i)elx M, 

for some constants C\, C 2 > 0. 

Let us introduce a triangulation T h of M, that is, a set of curved polyhedra K C M whose 
vertices are joined by geodesic lines. We assume that M — {J Ke r h K an d that f° r every distinct 
elements K\ , K 2 G T h , the set K\ n K 2 is either a common face of K\ , K 2 or else a submanifold 
with dimension at most [d — 2). To each element K we associate the set of its faces e by dK, and 
we denote by K e the unique element distinct from K sharing the face e with K. The outward 
unit normal to an element K at some point x G e is denoted by n ei K[x) G T X M . Finally, \K\ 
and |e| denote the d- and (d — l)-dimensional Hausdorff measures of K and e, respectively, and 
Pif = J2 e€dK |e| is the measure of the boundary of K. When N is a submanifold of M, we denote 
by dvN the volume element on N induced by the metric g. 

We impose that 

(5.3) h := sup d g (x,y) — ► 0, 

where d g (-,-) is the distance function associated to the metric 5. In particular, elements may 
become flat in the limit as h tends to zero. Interestingly enough, this condition on h is slightly 
more general than the corresponding condition in [2] , where the exterior diameter of the elements 
was considered. We denote the time increment of the discretization by r = t(/i), and set t n := rn 
for n = 0, 1, We also assume that, as h tends to zero, 

(5.4) t -> 0, /i 2 /t- -> 0. 

To each element if and each edge e G (9if we associate a family of locally Lipschitz continuous, 
numerical flux-functions f e> K '■ R x K — > K satisfying the following properties. 

• Consistency property : 

(5.5) f eK (u,u) = -i- [ g(f(u,x),n e . K (x)) dv e , ueR. 

|e| Je 

• Conservation property : 

(5.6) f e>K (u,v) = -f e ,K e (v,u), B,D6l. 
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• Monotonicity property : 

(5.7) lf e , K >0, ^f e , K <0. 

Note that, in the Euclidian case and when the flux depends only on the conservative variable, 
the consistency condition (|5.5p reads f e: K(u,u) — f(u) ■ n ej x, which is the usual consistency 
condition (c.f. [2]). For the sake of stability we impose the following CFL condition, for every 
element K E T: 

feJ<(u,v) - fe,K(v,v) 



TPK 

; —, max sup 

/g g\ |iv| K£T h ,e£dK u , veR ^ v 

T P K t • en / i 
-|^|- Lip(/) < 1, 

where 



< 1, 



Lip(/) := sup ^-(u,x) 

TIG? 



and | ■ \ g represents the norm induced on each tangent space T X M by the metric g. 

As in the Euclidian case, the finite volume method can be motivated by formally averaging the 
conservation law over an element K and using Stokes' theorem: 



= JtW\L u dVM + W\l K V9 ' fM dVM 

dvM + TFT / 9(f(u,x),n(x)) dv dK 



dt\K\J K — ' \K U9K 



~dl~\k\ J K udVM + \JC\ ^ Jj(f(%x),iL e>K (x))dv e . 



We then introduce the approximation 

U K~J]^J udv M, f e ,K{uK,U r Kj « j^j J g(f(u,x),U etK (x)) dv, 

and discretize the time derivative with a two-point scheme. 

Based on the above the finite volume approximations are defined by 

(5-9) u n K +1 :=u n K - £ T^Se,K{u n K,u^). 

The initial condition (|5.ip gives us 

(5.10) u° K := — { u (x)dv M - 

\ k \Jk 

Finally we define the function u h : R + x M -> R by 

(5.11) u h (t,x) = u n K , xeK, t n <t<t n+1 . 

We can now state our main convergence result. 

Theorem 5.1. Consider the Cauchy problem (|l.ip - (15.1|) associated with a conservation law on the 
Riemannian manifold (M,g) and some initial data uo € L l (M) n L°°{M). Suppose that the flux 
satisfied the growth condition (|5.2p . Let T h be a triangulation and t — r(h) be a time increment 
satisfying the conditions (|5 ,3|) - (|5 .4[) . and let f & ,K be a family of numerical flux-functions satisfying 
the consistency, conservation, and monotonicity conditions (|5.5p - (|5.8p . Let u h : M+ X M — ► M 
be the sequence of approximate solutions constructed by the finite volume method (|5.9p - (l5.1ip . 
Then, for all T > 0, the sequence u h is uniformly bounded in L°°([0, T), L 1 {M) n L°°(M)) and 
converges almost everywhere towards the unique entropy solution u £ i°°(R + ,L 1 (M) n L°°{M)) 
of the Cauchy problem ()l.ip - ()5.ip . 

The proof of this theorem is the subject of the next section. 
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Remark 5.2. Without loss of generality, we may assume that Lip(/) > 0. Indeed, in the trivial 
case Lip(/) = 0, f does not depend on u and the conservation law (|1.1[) has the explicit solution 
u(t,x) — uq(x) — £V g -f(x). Thus, for any triangulation, if we set f e K '■— j^j J e g{f{x),i\{x))dv e , 
we immediately find from (|5.9[) that u\ = u\ — nr'Vg-f(x) for all K. This gives u h (t,x) — 
u h (0, x) — tV g -f(x), and so lim^o u h (t, x) = uq(x) — tV g ■ f(x), which is the exact solution. 



6. Proof of convergence 

This section contains the proof of Theorem 15.11 We follow the strategy of proof proposed by 
Cockburn, Coquel, and LeFloch in the Euclidian case [2]. As we will see, the derivation of the 
required estimates for a conservation law on a manifold is significantly more involved. 

6.1. Discrete maximum principle, entropy dissipation and entropy inequalities. We 

begin with the discrete maximum principle and the L 1 contraction principle. 

Lemma 6.1. The finite volume approximations satisfy 

(6.1) max \u%-\ < ( max \u° K \ + dt n ) e° 2t " , 

and, given another approximation sequence v h (defined in the obvious way from a different initial 
data v E L l {M)^L°°{M)), 

(6.2) \\u h (t + r)-v h (t + T)\\ Ll{M) < \\u h {t)-v h {t)\\^ {M) . 

Proof. Let us first point out the following discrete version of Gronwall's inequality: if w n is a 
non-negative sequence satisfying Wq < a and 

n-l 

(6.3) w n < a + 6 / J Wj, n > 1, 

3=1 

for some non- negative constants a, b, then it follows that 

(6.4) w n <ae bn . 
Namely, let us define (p n by <fto — a and 



a + bJ2 

3=1 



n-l 



We have <j> n — <p n —i = &(/>„_ i, or 4> n = (1 + b)4> n -\. By induction we get cf> n = a(l + b) n , and the 
desired estimate now follows from 1 + b < e b and w n < cf> n . 

Let us turn to the proof of (|6.1[) . Using (|5.9|) and the consistency condition (|5.5p . we easily 
write 



„+l = ( 1 + J_ V- h,K(u\,u n K ) - f^ K {u n K ,u\) 
K \ \K\ ^ u n K -u n K 1 

-\k\ ^r^ K \^- Wl j K ^-f(u K ,x )dVM . 



e£dK K ? K 

The monotonicity condition (|5.7[) implies 



2-^ „n _ ,.n - ' 



so that, using the CFL condition 



w\ ^ ^t^k w p(/) - 
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Thus, we see that u 7 ^ 1 is a convex combination of u\ and (u 7 ^ ) e e9A, up to a divergence term 
in (16. 5p . Therefore, we can write 



% +1 < max(uK, max4 ) - t— , / V g -f(u] c ,x)dv M , 

eeoK XV 



A' 



^>min(<,min<J-^/ A ; 
and, by induction, 

n-l 1 
^ max I u< aI + E t max T17T / \^g-f( u K> x )\ dv M 



n-l 

< max \u° K \ + > t max (C\ + Coln-LI) 

3 =0 



n-l 

< max \u°fc\ + mC\ + > TC2 max lu^i 

KdT h *r~L AeT h ,v 



for all n > 0. Using GronwalPs inequality (|6.3[) - (|6.4p we arrive at the desired inequality (|6.1[) . 

We now turn to the proof of the L 1 contraction property (|6.2|) . It will follow immediately once 
we check that 

(6.6) Yl K +1 -< +1 H^I< E \u n K~vl\\K\. 

KeT h KeT h 

Introduce the discrete solution operator H, 

u n+1 ee {u n K +1 ) KeTh = {H K {u n )) KeTh =: H{u n ). 

Using the CFL and monotonicity conditions (|5.8p and (|5.7[) . it is easy to check that H is a 
monotone operator: 

(u n K <v n K , K e T h ) => (H K (u r x) < H K (v]\), KeT h ). 

Setting a A b = max(a, b) and a V b = min(a,6), we have Hk{u\\ V v]\) = V and 

H k (uk A v\\) = u]^ 1 A v^ +1 - Next, using the definition ([O]) we find 



r|e 
W 



if v "A' ~ "A /x y A - "A v U K - a K ^ U K 

E W\ (f^( U K V <> « V "SW - A V K' « A <)'))■ 



e£dK 1 

Observe now that 

Ka V w£) e < u n Ke V < c and « A v n K ) e > u n Kc A v\ c 
This, together with <9/ e ,A/<9v < 0, implies that 

/ e ,A« V V%, « V t£) e ) - / e ,A« A « A V%)e) 

> fe,K(u n K V V n K ,U n Ke VVWJ- fe.K(u n K A A t&J, 

so that 

U A +1 V W A +1 - U A +1 A W A +1 < u ^ V w " A 



E iS V <, < e V t&J - U« A „£, u n Ke A 



eGOA 1 



The inequality (|6.6|) now follows by multiplying by |-ftr| and summing in K G T , since the 
conservation property (|5.6p allows us to cancel all of the numerical flux terms. The proof of 
Lemma |6. II is complete. □ 
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The discrete entropy inequalities will be formulated in terms of the following convex decompo- 
sition of u\. For each K, e, define (cf. |2J, formula (3.6)) 

(6-7) U n + e ' := U\ - ^(f ejK (u K ,U K J - feM^ul)) 

and set 

(6-8) := u n K +1 e -Tk\I k V ^«' x) dv M . 

In view of (|5.9p and the consistency property (|5.5p . we have 

(6.9) z4 +1 = ± £ |e| tt n+i. 

PK el^K 

The following lemma establishes the existence of numerical flux terms and a first entropy 
inequality relating and u K . We omit the proof and refer to [2] and the references therein. 

Lemma 6.2. Let (U,F) = (U(u),F(u,x)) fee a convex entropy pair. Then there exists a family 
of numerical flux function F e> K '■ R 2 - * R satisfying the following conditions: 

• F ei K is consistent with the entropy flux F: 

F etK (u,u) = -j-r [ g(F(u,x),n e . K (x)) dv e , M6l 

M J e 

• Conservation property : 

Fe tK (u,v) = -F e! K e (v,u), li,«6R. 

• Discrete entropy inequality: 

(6.10) u{u n +l) - c/«) + ^(v(4,4 e ) < o. 

We can now write the discrete entropy inequality in terms of it'^v and u^- (see (I6.10[) ) : 

(6.11) l^l V ' 

The following entropy dissipation estimate is the key to the proof of convergence, as was 
observed in the Euclidian case by Cockburn, Coquel, and LeFloch [SJ. 

Proposition 6.3. Let U : R — > R fee a strictly convex, smooth function. Then, for all n we have 

E + f E ™Ke x -< +1 l 2 

( 6 - 12 ^ 1,11*1 

< E E r\e\F e>K (u n K ,u n K )+ E ^ U ^ e 1 , 

KgT h /fgT k KeT h PK 

eedK eedK eedK 

where a denotes the modulus of convexity ofU. 

Proof. First, we observe that from an elementary lemma on convex functions (Lemma 3.5 of [2j) 
and (|6.9[) one easily obtains 

(6.i3) e WK +1 ) + f E ^\<l-< +l \ 2 < E ^u{u^). 

T ~L, * t —\. PK , PK 

KeT h KeT h KeT h 

eedK eedK 

Next, we sum up the inequality (|6 . 1 1 [) with respect to all K and e, and divide by px to obtain 
(6-14) E — UVgl)- E U(ul)\K\- E T\e\F e , K (u K ,u K )- E < 0. 

— PK — — — PA" 

Arer' 1 KeT h KeT h KeT h 

eedK eedK eedK 
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Note that, again, all of the flux terms F e ^ {u\, u\ ) cancel out due to conservation. Combining 
(|6.14p with (|6.f 3[) gives (|6.12|) . The proof of Proposition 16. 31 is complete. □ 

As a consequence of the entropy dissipation estimate (|6.12p , one can deduce the following 
bound: 

(6-15) Yl E ™Ke-<| 2 <\\uo\\h {M) +C(T), 

n=1 KeT h PK 
eedK 

where C(T) is a constant depending only on T and N is an integer such that T — Nt < r. For 
instance, it suffices to take U(u) — u 2 /2 in (|6.12[) . 

We now derive a local entropy inequality which does not involve the numerical flux and is a 
consequence of the inequality (|6.1ip . The proof is completely similar to that of Lemma 3.6 in [2] 
and we omit it. This inequality will then allow us to prove a global discrete entropy inequality. 

Lemma 6.4. Let (U, F) be a convex entropy pair. Then for all elements K and all faces e £ dK , 
we have 

/g_ 16 J PK PK PK? PK? 



+ f [ g(F(u'f C ,x) — F(uf C) x),n e K(x)) dv e < ^R n +l + —R^l 

\ e \ Je PK ' PK e 



where R^ e is defined in (|6. 1 1|) . 

Proposition 6.5. Let (U,F) be a convex entropy pair and let </> = 4>(t,x) £ C c ([0,T) x MJ be 
test function. For each element K and each face e of K , set 

(6.17) 4>™ = -j-r f " + 7 0(t,x)dv e dt, 4> n K =y2 -^c 

r\e\Jt n Je e edK PK 

(6.18) §^ n = -(4> n K -4> n K - 1 ). 

T 

Then, u\ satisfies the following inequality: 



(6.19) »-"A,r" 



EE f " +1 f (u{u n K )d^K n +g(F(ul,x),V g cj > (t,x)))dvMdt 

n K(7 h " 

+ E / U(u° K )4>° K dv M <E h + Q h + R h , 

JK 



KeT h 

where the error terms E h , Q h and R h are given by 

(6.20) £ k = EE ^^KTeX^r-C), 

eedK 



N 

Q h -EE Qk, 

n=0KeT h 



(6.21) 



Qk-=Y. t 1 9{F(u n K ,x),n e , K (x))dv e ^- f" 1 f g{F{u n K ,x),V g( t>(t 1 x)) 



eedK 



dvxdt, 

t n JK 



& = E E ™«*s# 

(6.22) r l= o KeT » PK 

eedK 

Rte=U(u n +t)-U(ul + t). 
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Proof. Multiply the inequality (|6 . by 4>™\e\ and sum over all K 6 T , e s Note that the 
first four terms in (|6.16[) . as well as the last, give the same result whether they involve K or K e 
(since the sum is taken over the same set of elements). We thus obtain 



2 E 

KeT h 
eedK 

+ E T MC(n [9(F(u n KB ,x)-F(u n K ,x),n e<K (x))dv e )<2 ]T ^^C- 
Ker , Me|J e Kerfe PK 

eedK eedK 

Also, due to the conservation and consistency properties of F et K we have 

E FT [ 9(F(u'J c ,x),n e>K (x)) dv e = - V i f g(F{u% ,x),n e>K (x)) dv e , 



which implies 

ifGT h ^ 

(6.23) eedK 

E r N^n fg{F(u K ,x),n e , K (x))dv e < — *0?. 

KeT h l6lJe KeT h PK 

eedK eedK 

Observe now that, setting Q K as in (|6.21|) . we deduce from the definition of 0™ (see (|6.17|) ) that 

Ker* |e| L 

(6.24) ee9K 

= E /" 1 / ff(J , («x»*)»V0(t,a;))d» w dt+ J2 ' 

KeT h Jt n JK KeTh 

From the convexity of U and Jensen's inequality, we see that 

(6.25) \K\U(u n K +1 Wk < E T^^Ofe- 

eedif 

Thus, from (|6.23[) . (|6.24p and (|6.25[) and after summation in n, we obtain 

J2 E \K\(U(u n K +1 )-U(u n K ))$ n K -Y, E / / 9(F(u K ,x),V<f>(t,x))dv M dt 

n=0 KeT h n— KeT h K 

n n , 

<^ + E E « + E E 

n=0KeT h n =° KeT h 

eedK 

with £ given by (|6.20p . Summing by parts in the first term, one obtains the desired entropy 
inequality (|6.19[) . The proof of Proposition 16.51 is complete. □ 

6.2. Measure-valued solutions and convergence analysis. We now prove the strong con- 
vergence of the approximate solutions u h given by the finite volume method towards the entropy 
solution to the Cauchy problem under consideration. To this end, we need to rely on the frame- 
work of measure-valued solutions to conservation laws developed in [4] in the Euclidian case and 
extended to manifolds by Ben-Artzi and LeFloch pp. The key question is to justify the passage 
to the limit as h — > in the inequality (|6.19|) and use the entropy dissipation estimate ()6.15j) to 
control the error terms. 

It is well-known that to the sequence u h (which is uniformly bounded in L°°([0,T) x M) for 
every T > 0) we can associate a subsequence and a Young measure v : [0, T) x M — > Prob(R), 
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which is a family of probability measures in R parametrized by (t,x) G [0, T) x M. The Young 
measure allows us to determine all weak-* limits of composite functions a(u h ), for arbitrary real 
continuous functions a, according to the following property : 

(6.26) a(u h ) (v, a) when h -> 

where we use the notation 

(v, a):— / a(A)dz/(A). 



In view of the above property, the passage to the limit in the left hand side of (|6.19p is (almost) 
immediate. The key uniqueness theorem [3J [T] tells us that once we know that v is a measure- 
valued solution to the conservation law, we can prove that the support of each probability measure 
Vt,x actually reduces to a single point u(t, x), if the same is true for t = 0, that is, v t . x is the Dirac 
measure S u n )X y It is then standard to deduce that the convergence in (|6.26|) is actually strong, 
and that, in particular, u h converges strongly to u which in turn is the unique entropy solution 
of the Cauchy problem under consideration. 

Lemma 6.6. Let v : [0, T) x M — > Prob(R) be a Young measure associated with the sequence 
u h . Then, {or every convex entropy / entropy-flux pair (U, F) and for every non-negative, smooth 
function 4> : [0, T) x M — > with compact support, one has 
(6.27) 

-/ / (u t , x ,U{-))d t (j)(t,x)+g((ut, x ,F(;x)),V(l)(t,x))dvMdt+ [ U(u (x))<t>(Q,x) dv M 

JO JM JM 

< [ I (vt, x ,V g -F(;x))(t>(t,x)dvMdt- [ [ {Vt,x,U'(-)V g -f{;X))(t>{t,x)dVMdt. 



JO JM JO JM 

Proof. First, we claim that the discretizations (|6.17p . (|6.18|) satisfy the following estimates, whose 
proofs are elementary: 

sup \d t <t>{t,x) - d t (f>k | < (r + h K )\\4>\\ C 2 {[tnitn+l]xK) , 
te[t n ,t n+1 ] 

and, for each face e, 

\€ - 4>k\ <(t + h K )U\\c^[t n ,t n+1 ]xK)- 
In view of these estimates and (|6.26l) , the left hand side of (|6.19[) converges to the left-hand side 

of Km . 

We now consider the error terms on the right-hand side of (|6.19l) . First, note that E h equals 
the corresponding one in [5] , so we refer the reader to that paper for a proof that E h converges 
to zero. Note that the use of estimate (|6 . 1 5[) is essential in that proof. 

Next, let us prove that 

(6.28) \imR h = -[ [ {ut, x ,U'(-)V g -f(;x))<l>(t,x)dvMdt. 



h— >0 



J M 



Recall that R h is given by (I6.22|) . Introducing the averages 

4>(t, x) dvudt, 



I rtn+l 



VK - W\r J, 
we claim that it suffices to prove that 



K 



(6.29) lim (R h U '( u k) t f V 9'/«: =0 dv M <t> n K ) = 0. 

n KeT h K 
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Indeed, we have 



E E u 'Wk)r I y g -f(u n Kl x)dv M r K 

n K£T h K 

= -EE U 'Wk) / ^ g -f(nl,x)cf > (t,x)dv M dt 

„ TS^rrl, Jt„ JK 



n KeT h 



E E U '( U K)(j t n+1 j K V 9 -fWk,x)dv M ^, -J <Kt,x)dv M dt 



n K<ET h 



lt„ JK \^\JK 
rtn+i 

/ ^g-f{u n K ,x)4>{t,x)dv M dt), 



so that taking the limit as h — > we find, for the first term, 



lim- / / U l (u h )\7 g -f(u h ,x)(j)(t,x)dvMdt 

lo JM 

[ fax, U'{-) V g -f(;X))(f>(t, x) dv M dt, 
lO JM 

while the second term easily converges to zero. 

We now prove (|6.29[) . First, observe that from the definition of </>^, it is easy to see that we 
may as well consider the term R h with 0^- instead of 0™. We thus have 

N 

Vg-f(uK,x)dv M 4>K 



9 
IK 



'JM 



eedK 

eedK 

+ ™C/'«)^ / V 9 -f(u K ,x)dv M( f ) % 

eedK 

= E E ™^(^(<e 1 )-^e 1 ) + ^K +1 )^/ V s ./K +1 
n=0 PK ' ' 



i KeT" 

eedK 



+ E E ^ f/ '(^ +1 )^ / v 9 -/(^ +i ,^^m(^ +1 -^) 



™=0ifeT h 

eGdif 



- E U '( u °k)t I V g -f(u° K ,x)dv M <f>K = A + B + C. 

Note that all sums in n may be written from to N because of the compact support of <p. 

Consider the first term, A. From the bound (|5 . 2[) on the divergence of / and from the maximum 
principle (|6.1j) we conclude that 

^ r Jv g -f{u n I t\x)dv M = 0{T). 
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Thus, Taylor developing U around the point u^, 1 , and using (|6.8p , the definition of <+\ we find 

I^E E ^|f/ , (< e 1 )0(r) + 0(r 2 )-f/'(< +1 )0(r)|^ 
n xer h 

Next, using the estimate f|6. 1 5[) and the Cauchy-Schwartz inequality, we find that the term inside 
the sum above is bounded, so that A = 0(t). 

Consider now the second term, B. The function <j) is regular, so that </>^- +1 — 4>\ = 0(t). 
Therefore, we have 

1*1 ^E E \K\0(t*) = 0(t). 

n K£T h 

Finally, using again the bound (|5.2p , we see that the last term C tends to zero with t. This 
concludes the proof of (|6.28[) . 

Next, consider the term Q h = J2 n k Qk gi ven by (|6.2ip . Using Green's formula, we can write 



Q'k = E Tl / 9{ F ( u K, x ), n e : K{x)) dv e [ [ <j){t, x) dv e dt 

^ ai7 l e l Je Jt„ Je 



eedK 



g{F{u K ,x)(j){t, x),n(x)) dv aK dt + I f V g -F(u r ^,x)(f>(t 7 x)dvMdt. 

dK Jt n JK 

We wish to obtain 

(6.30) limQ h =[ [ {v t , x ,V g -F(;x))(f>(t,x)dvMdt, 

Jo Jm 

so that according to the above equality we must prove that 

lim ^2 (|~| J 9{F{u%,x),n etK {x)) dv e £ J^cf>{t,x) dv e dt 



n KeT h 
(6.31) eidK 



g(F(u K ,x)4>(t,x),n e , K (x)) dv dK dt) = 0. 



Let x e be any point on the edge e. Note that 

rh Jt„ 



n KeT 
eedK 



E E J J^g{F{v?K,x),Ti etK {x)) Q- Jj>(t,x)dv e - cf>(t,x)^jdv e dt 

dK 

E E J t J^(9{F{u K ,x),n e ,K(x)) -g(F(u K ,x e ),n etK (x e ))y 

■ (A I <f>(t,x)dv e - 4>(t,x)^Jdv e dt, 



KeT n 

eedK 



so that taking absolute values and in view of the regularity of F and <f>, we can bound the last 
term by 

CT J2 VKh 2 . 

KeT h 
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using the CFL condition (|5.8[) . we can finally write 



Q h - ( ( V g -F{u h ,x)<j>{t,x)dv M dt 
Jo Jm 



^ w ' KeT h 

< c y = —c\m\, 



T T 

KGT h 



so that taking the limit as ft. — > and using (|5.4|) yields 



T 



Q h -> / / (v t ,x,V g -F(;x))(f>(t,x)dv M dt, 



M 

which is (|6.30p . This completes the proof of Lemma [6T6J □ 

Remark 6.7. From the proof of the previous lemma, we see that if the divergence-free condition 
(|1.2[) is satisfied, then the error term Q h converges to zero with h and R h vanishes altogether. In 
that case, we obtain, instead of (|6.27p . the usual weak formulation 

f (v t , Xl U{-))d t 4{t,x)+g((vt,x,F(-,x)),\7(f)(t,x))dvMdt+ f U(u (x))^(0, x) dv M < 0. 
Jm Jm 

Proof of Theorem \5.1\ According to the inequality (|6.27p , we have for all convex entropy pairs 
(U,F), 

d t {v,U{-)) + V g -{{v,F{.,x))) < (v, (V g -F)(-,x)) - (v,U'(-)(V g -f)(;x)) 
in the sense of distributions in [0, T) x M. Since, for t = 0, the Young measure v is the Dirac 
mass S Uo (because uq is a bounded function), we know from [T] that there exists a unique function 
u € L°°([0, T) x M) such that the measure v remains the Dirac mass 6 U for all time < t < T . 
Moreover, this implies that the approximations u h converge strongly to u on compact sets at least. 
This concludes the proof. □ 
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